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Abstract 

We examine how nonperturbative effects in string theory are transformed under 
Q_i! the T-duality in its nonperturbative framework by analyzing the c = 1/2 noncritical 

string theory as a simplest example. We show that in the T-dual theory they also 
take the form of exp(— So/g s ) in the leading order and that the instanton actions So 
of the dual ZZ-branes are exactly the same as those in the original c = 1/2 string 
^ ■ theory. Furthermore we present formulas for coefficients of exp(— So/g s ) in the dual 

1 theory. 



1 tkuroki@post.kek.jp 

2 fumihiko _sugino@pref . okayama. j p 



1 Introduction 



T-duality is a relation between perturbative vacua in string theory and is known to hold at 
each order in perturbation theory for the critical string theory pQ. Since it is a character- 
istic feature of string theory, it is believed to play an important role even in constructing 
its nonperturbative formulation. From this point of view, it is quite intriguing to ex- 
amine how nonperturbative effects are transformed under the T-duality in a possible 
nonperturbative framework of string theory. As such a framework, the noncritical string 
theory provides a useful toy model, because it can be formulated nonperturbatively by 
matrix models and their nonperturbative effects are identified as a kind of instanton ef- 
fects [2]~[Z]- Moreover, we can explicitly formulate the T-duality in the noncritical string 
theory, for example, in the c = 1/2 string theory as discussed in [8]. 

However, as for the c = 1/2 noncritical string theory, which does not have a continuous 
target space, it has been shown explicitly in [8] that the theory is not invariant under 
the T-duality for lack of global winding modes associated to string world sheets of higher 
genus topology. This viewpoint for the T-duality still stays at the perturbative world sheet 
picture, although the corrresponding string field theory has been explicitly constructed. 
Therefore, as a first step to understand the T-duality truely at the nonperturbative level, 
it is important and interesting to identify nonperturbative effects in the dual c = 1/2 
string theory and to clarify how they correspond to those in the original c = 1/2 string 
theory. This is the aim of this paper. 

Since as shown in refs. [HIE] the dual theory is defined by the 0(n) model on a random 
surface [TU] with n = 1, we analyze this model and show that the leading nonperturbative 
effects in the dual theory also take the same form of exp(— So/g s ) as in the original 
theory, where g s is the string coupling constant. This form of the nonperturbative effects 
in general reflects the large-order behavior of perturbation series in string theory [TT] . 
Thus, we find that the perturbation series in the dual theory also shows the same large- 
order behavior. Moreover, for the standard c < 1 noncritical string theory, the values 
of So are deduced from the string equations [21 [T21 EG] , or directly from matrix models 
as their instanton effects mentioned above. Other techniques deriving So are also found 
in refs. [T¥| [15]. For the 0(n) model with general n on a random surface, in contrast 
to the matrix models above, the string equation has not been derived (to the best of 
our knowledge), from which the nonperturbative effects in the dual theory can be seen. 
However, since the total free energies of the original theory and the dual theory are 
equivalent by definition for the case n = lH, the free energy of the dual theory should 
satisfy the same string equation as that of the original theory does. Therefore, we expect 
that the dual theory has exactly the same values of So as those of the original theory. 
We will see that this is indeed the case. Since it is known in the standard noncritical 
string theories that So can be identified [16]~[20] as the classical actions of the ZZ-branes 
|21j . we can conclude that the classical actions of the dual ZZ-branes are the same as 
those of the original theory, and it gives an important basis in determining the T-duality 

3 Because of the subtlety of wcll-dcfmedness for the matrix integrals in the double scaling limit, we 
can safely say that the equivalence holds in the sense of the 1/iV-expansion. 
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transformation rule of the ZZ-branes. 

The organization of this paper is as follows. In the next section, we formulate the 
T-duality in the c = 1/2 noncritical string theory in terms of the two-matrix model, 
and we see that the dual theory is defined by the 0(n) model on a random surface with 
n — 1. In Section 3, we make the saddle point analysis in the 0(1) model and derive 
nonperturbative effects in the dual c = 1/2 string theory. We proceed to examine the next- 
to-leading contribution to these nonperturbative effects in Section 4. The last section is 
devoted to the conclusions and discussions. In Appendix A, we explain the double scaling 
limits of the original and dual two-matrix models, and present a derivation of resolvents 
in the dual model. Some other computational details are discussed in Appendix B. 



2 T-duality in the c = 1/2 noncritical string theory 

The original c = 1/2 string theory is defined by the double scaling limit [22] of a two- 
matrix model 

Z = j d N2 Ad N2 Bexp[-NtrS(A,B)}, 

S(A,B) = -A 2 + -B 2 — cAB — -A 3 — -B 3 , (2.1) 
v ' ; 2 2 3 3' V ; 

where A, B are the N x N Hermitian matrices. This model describes the Ising model on 
a random surface [23] with the inverse temperature given by 

P = ~]ogc (2.2) 

Various amplitudes in the c = 1/2 string theory can be computed based on ( 12. ip . for 
example by solving the Schwinger-Dyson equations as discussed in refs. [24"1 125] . 

In order to perform the T-duality transformation to this model, we first note that 
the familiar T-duality transformation in closed string theories compactified on a circle is 
basically the same as the Kramers- Wannier dual transformation to the XY-model defined 
on the string world sheet (see, e.g. [26] or Appendix A in [8]). From this point of view, 
the T-duality transformation for this model will be naturally formulated as the Kramers- 
Wannier transformation of the Ising model on a random surface. Since the matrices A 
and B can be regarded as the up and down spins on a random surface, respectively, the 
Boltzmann factors of the original Ising model on a random surface are related to the bare 
propagators for A and B as 

( A v A ki) haxe = ( B v B ki) haxe = -^SuSjkLe* 3 , (AijBu)^ = ^5 u 5 jk Le-^ (2.3) 

with L = \fcj{\ — c 2 ). Therefore, the Boltzmann factors in the dual model should be 
proportional to e ±/3 , where j3 is obtained by the Z 2 Fourier transformation 

= K(eH + e-P), 
e -P = K(eP-e-P). (2.4) 
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K is an overall normalization constant given by K = (e 2/3 — e 2 Py/ 2 . It is easy to introduce 
matrices realizing such bare propagators. If we define new matrix variables 

X = -L(A + B), Y = -L(A-B), (2.5) 

then their bare propagators exactly give the Boltzmann factors of the dual model: 

(X i5 X kl ) haxe = Ls a 6 jk -j==et, (YijY kl ) haie = — <5 a <5 ife -^J==e^ '. (2.6) 

In this sense, these represent "stick" or "flip" of the dual spin, respectively. Thus, the T- 
duality transformation amounts to making a field redefinition ( 12.5ft and the dual c = 1/2 
string theory is defined by the double scaling limit of a dual two-matrix model [HI E] 



Z = j d N2 Xd N2 Y exp[-Ntr S D (X,Y)}, 



S D (X,Y) = i^X 2 + i±^F 2 - | (X 3 + 3XF 2 ) , (2.7) 

where g = gj y/2. This model is also known as the n = 1 case of the 0(n) loop gas model 
on a random surface |10j . 

The T-duality transformation (12.51) is a trivial change of the integration variables, and 
the total free energies given by 

* r = -Jp*> S Z, (2.8) 

take the same value at least order-by-order in the 1/X-expansions for both models. How- 
ever, the coincidence of the free energies does not always lead to the T-duality. Indeed, we 
can see that the T-dual relation is broken in correlation functions on a surface of higher 
genus topology. For example, it is shown explicitly in [S] that disk amplitudes in both 
theories have the same functional form, while a disk amplitude with one handle in the 
dual theory have different functional form from a corresponding amplitude in the original 
theory. More precisely, for the universal part of the disk and cylinder amplitudes, the 
following identification between the original and dual models holds: 



7l( tr <V tr <^) 



tr^, (2.9) 



in the sense that they take the same functional form as functions of ( and £ respectively 
with certain identifications of parameters. This is expected because both operators in 
(I2.9p are interpreted as the Dirichlet type boundary conditions for the original and dual 
spins respectively, under the T-dual relation. Recall that (XjjXfc/) bare represents "stick" 
of the dual spin. On the other hand, the disk amplitude with one handle 
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no longer has the same form as 




1 1 



> 



i 



(2.11) 



even for the universal part. (The subscript "1" put to the expectation values represents 
the random surface having one handle on which the expectation values are evaluated.) 
This difference originates from the excitations of odd number of F-loops along a handle. 
Note that, if they are present, such a configuration cannot be interpreted as a dual spin 
configuration because (^•^) bare represents the flip of the dual spin. Of course, from 
eqs. (12.51) we have 



However, it does not imply the T-dual relation between the amplitudes, because the 
operator appearing in the r.h.s. is not consistent with the interpretation as "stick" of the 
dual spin in the l.h.s. Thus, we should remark that, in order to show the T-duality, it 
is necessary to give a consistent interpretation of X and Y as dual spins, not only the 
transformation of variables (12. 5p . It is the same situation to the case of the Ising model 
on the regular lattice. The high temperature expansion of the Ising partition function 
on the plane has an one-to-one correspondence to the low temperature expansion. (For 
example, see [27].) It is a manifestation of the Kramers- Wannier duality (the T-duality). 
All the terms in the high temperature expansion can be written as configurations of loop 
gas. In the case of the surface with higher genus, some of them contain the loop gas 
configurations where odd number of the loops surround topologically nontrivial cycles. 
They can not be interpreted as dual spin configurations and violate the T-dual relation. 
Therefore, we can conclude that, although the total free energies are same between the 
original and dual theories, the T-duality does not hold in their higher genus parts due to 
the excitations of odd number of y-loops along the handles. 

In string theory with continuous target space, the T-dual symmetry arises when the 
target space is compactified. It is a symmetry under the interchange between momentum 
modes and winding modes in the compactified directions. From this viewpoint, the Ising 
model, whose target space is discrete (consists of two points), contains counterparts of 
the momentum modes, but not those of the winding modes for either case of random 
or regular lattice. We can understand that this asymmetry between the momentum and 
winding modes is the origin of the breaking of the T-dualityLl 

In the original model given by ( 12. ip . we can apply the method of the orthogonal 
polynomial [29J and derive the string equation in the double scaling limit as 




4 It is possible to introduce the counterparts of the winding modes to matrix models with discrete 
target space to have the exact T-dual symmetry as discussed in [25] , 




(2.12) 




(2.13) 
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where t is the cosmological constant, g s is the string coupling constant, and /(£) is the 
second derivative of the free energy fit) = g 2 s F[t) as a function of t [301 l3"lffi From this 
equation we can deduce the asymptotic expansion of F as 



9 t* 1 , 
FW = 28^J + 24 h * t + "' 



(2.14) 



which is nothing but the genus expansion. Since as shown in Appendix |A] (eqs. ( 1A.4I) . 
(1A.6P ) the double scaling limit is taken for the dual model (12.71) in the same way as in 
the original model [8], the free energy of the dual model should also satisfy the same 
string equation as in (12.131) in the double scaling limit@ and should have the same genus 
expansion as given in (12.141) . However, even in this case, it is possible that they have 
different nonperturbative effects. Namely, suppose fi and f 2 satisfy the same string 
equation (12.131) . then the semi-classical treatment of (12.131) with g s ~ h leads to the 
difference Af — fi — f 2 of the form: 



which yields a nonperturbative ambiguity of the free energy F(t) in the following form: 



C\ and C 2 are numerical constants which cannot be determined from the string equation 
alone. Therefore, if we interpret one of j\ and f 2 as a quantity of the original model 
and the other as that of the dual model, as long as C\ or C 2 is nonzero, the dual theory 
has a different nonperturbative effect from that in the original theory. It is known that 
the nonperturbative effects in the original theory itself take the form ( 12.161) [31 [7] , so the 
dual theory must also have nonperturbative effects of the same form (up to the overall 
constants). We will see that this is indeed the case by computing the nonperturbative 
effects in the dual theory directly from the matrix model (12. 7p . It is well known that the 
exponents of the nonperturbative effects are identified as the actions of the ZZ-branes and 
are provided by the disk amplitudes in the presence of them. Therefore, the fact that the 
dual theory also has the nonperturbative effects as in (12.161) implies that the actions of the 
ZZ-branes in the dual theory (the dual ZZ-branes) are the same as those in the original 
theory. It is worth noticing that the string equation can fix not only the exponents in the 
nonperturbative effects but the power of t in the factors in front of them. 

5 Note that we have changed the sign of the free energy compared to that in [3D] and multiplied it by 
2 because the potential in (|2.1|) is not even. 

6 In the dual model (|2.7|) . we have the expression of the integrals over the eigenvalues of X and Y 
(|3.4|) after integrating out the angle variables. It contains the factor Yii<j VCm* + MjOj m addition to the 
Vandermonde determinants usually appearing in the case of the original model. Due to the additional 
factor, the orthogonal polynomial method does not work well, and it makes directly deriving the string 
equation difficult. 




(2.15) 




(2.16) 
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On the other hand, the coefficients C\ and C2 in the nonperturbative effects cannot 
be fixed from the string equation. However, it is shown in [H [7] that, if we compute 
these coefficients in the c < 1 noncritical string theory directly from a matrix model, they 
turn out to be unique and universal in the sense that they do not depend on details of a 
potential in the matrix model. Thus, it is interesting to examine whether this is also the 
case with the nonperturbative effects in the dual c = 1/2 string theory and whether they 
agree with those in the original theory. Moreover, since the nonperturbative effects are 
computed from certain disk and cylinder amplitudes [H E] , it will be possible to find out 
the T-duality at the nonperturbative level from the knowledge of the T-dual relation of 
the disk and cylinder amplitudes. It is expected that the analysis reveals the existence of 
large universality including T-duality for nonperturbative effects in string theory, or for 
string theory itself. 



3 Nonperturbative effects in the dual theory 

In this section, we derive the leading part of the nonperturbative effects in the dual 
c = 1/2 noncritical string theory directly by evaluating instanton contributions in the 
matrix model (12.71) . 



3.1 Chemical potential of instanton 

As a preparation for instanton calculus in the dilute gas approximation, here we formulate 
the chemical potential of an instanton in the dual two-matrix model (12. 7h . 
By rescaling the matrices, the partition function in (12.71) becomes 



Z N {h) 
S D (X,Y) 



dXdY exp 
1 -c 



N 

--trS D (X,Y) 
n 



-X 



1 , \±^_y2 _ 1 , X 3 + 3X y2X 

2 3 K ' 



(3-1) 



where h = g 2 and the measures dX, dY are defined with the normalization: 



dX exp 
dY exp 



N 
~h 



tr 



1 -c 



N 



-X' 



1 + C y2 



(3.2) 



so that the free energy has the standard 1/iV-expansion 

F = -jp log Z N (h) = F (h) + ^F x {h) + ±F a (h) + 



(3.3) 



where Fk(h) represents the contribution from random surfaces with fc-handles. 
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The formula given in [2H] enables us to rewrite the partition function in terms of the 
eigenvalues of X and Y as 

z N {h) = d n ^ f (t[ d\d^ A n ^ A +^f ex P -f E + - A ^ 2 ) 

(3.4) 

where Aj, \ii (i = 1, • • • N) are the eigenvalues of X and K respectively, and A^(A) = 
YlN>i>j>i(^i ~ -\?) * s the Vandermonde determinant. V"(A) is defined as 

V(X) = i^A 2 - 1a 3 . (3.5) 

The proportional constant is explicitly computed in Appendix [B] (See eq. ( IB. 311) ). 

Hereafter, as a configuration of one instanton, we consider a situation where one pair of 
the eigenvalues (say (A^v, Hn)) is separated from the other pairs (Aj, /ij) {i = 1, • • • , N— 1) 
as a point on the (A, //)-plane. Setting x = A^, y = fiN, the partition function is expressed 
as 

J v \f = \ J rii<i<jv-i(^ + ^) il (y + A*0 

xe -f Ef = 7 1 (v(A ! )+i^M?-A 1 ^)-f(^)+i^^-V ) 

= ^ /^^_J^y- det( "T x/)de 5r y/) e -f»^V^> 
y y det(y + r / ) 

= -p^^-iC*) J d*dy{ ; e tt^( *) 

Z N -i(ti) J dxdy e~ v ^ x ' y) , (3.6) 



where X', V are (iV — 1) x (AT — 1) Hermitian matrices, and 



Z N ^{h!) = j dX'dY'e-^ s ^ x '< Y '\ 

(OY = - [ dX'dY'Oe-^ trS ^ x '> Y '\ (3.7) 

ZN-i\h') J 



In the above equation, we have defined h! as 

N N-l 



h " h! ' 

and V e g(x, y) as 



(3.8) 



-v eBix , y) = / det(x-X')det(y-Y') \'^ SnM 



det(y + Y') 



e -^> D (x,y)_ (3.9) 



As in the one-matrix model case considered in [J], the total partition function is divided 
into multi-instanton sectors as 

Z N (h) = Z { °-' mst \h) + Z { ^' mst \h) + • • • , (3.10) 

where the /c-instanton sector is characterized by fc-pairs of the eigenvalues separated from 
the other pairs on the (A, /i)-plane. More precisely, the partition functions in the 1- 
instanton sector and 0-instanton sector are given by 

Z^ iwt \h) = N^Z N ^(h') [ dxdye- v ^\ (3.11) 
Zt mst \h) = ^±Z N ^{h!) ! dxdye- v ^y\ (3.12) 



Dn J(x,y)eS 

respectively, where the factor N in (13.111) means the number of ways to specify a separated 
pair of the eigenvalues, and S is the support of the eigenvalue distribution on the (A, /i)- 
plane in the large- N limit. Taking the ratio between them, we obtain 

z (i-iMrt)/,x f dxdye- v ^ x ^ 
fi = Z f n . W =N f»W JL - (3.13) 

Z^ mst \h) f^ y)£S dxdye-y^y) 

Following the argument given in [I], it is easy to see that fi defined in this equation is in 
fact the chemical potential of the instanton, namely a statistical weight of the instanton, in 
the dilute gas approximation for the computation of the free energy of the dual two-matrix 
model. 



3.2 Saddle point analysis 



In this subsection, we apply the saddle point method to the numerator in f)3.13p . which 
is valid in the large- N limit. Since V e e(x,y) can be rewritten as 



-V e g(x,y) 



e tr log(x-X')+tr log(2/-Y')-tr ^g(y+Y') Y g - ^-S D (x,y) 



(3.14) 



one may expect that V e s(x, y) in the large- N limit can be expanded in terms of connected 
Green functions as 



-Veg(x,y) 



exp 



N-l 



h7 ^S D (x,y) + (trlog(x - X'))' d + (trlog(y - Y'))' d - (trlogfo + Y'))' 

+ \ <(trlog(x - X')f)' c + \ <(trlog(y - Y')) 2 )' c + \ ((trlog(y + Y')f)' c 
+ (tr log(z - X')ti log(y - Y% - (tr log(y - Y')tr log(y + Y% 



(tr log(x - X')tr log(y + Y% + 0[- 



(3.15) 
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where the subscripts "cf ' and "c" represent amplitudes of the disk and cylinder topologies, 
respectively. However, if x or y is inside the support of the eigenvalue distribution of X' 
or Y' respectively, the operator tr log(a; — X') or tr log(y — Y') becomes large by itself to 
make the expansion (13.151) not valid [32]. This motivates us to divide the numerator in 
(I3TT3D as 



X / dxdy e~ Vcs{x ' y) 

<(x,y)<tS 



X I dxdye~ v ^ y) 

dxdye~ v ^ y) + N j dxdy e~ v ^ x ' y) , 

(3.16) 



+X 



where S — Ix X Iy, an d Ix (W) is the support of the eigenvalue distribution of X' {Y'). 
In our choice of the potential, Ix and ly should be connected intervals, and their explicit 
forms are given in Appendix IA~1 as eqs. ( 1A.42I) and ( 1A.44I) . 

3.2.1 The first term in (IBTTBj) 

First we consider the case where both x and y are outside the supports of the eigenvalue 
distributions of X' and Y' . Then, the expansion (I3.15P is justified and the leading part 
of V e g(x,y) in the large- A" limit, denoted as v}^\x, y), is given by 



cxp 



- ^j^S D (x,y) + (trlog(z - X% + (trlogfo - Y% - (trlogfo + Y% 



exp 



N - 1 

-—^S D (x,y) + (trlog(x-X% 



(3.17) 



where in the last line we have used the Z 2 symmetry under Y — > —Y of the action (13. 1|) . 
Therefore, we have 



V^{x,y) = ^[V{x 



(v(x) + ^y 2 -xy^-(tT\og(x-X%. (3.18) 



Since V^\x, y) is proportional to N — 1, we can apply the saddle point method to evaluate 
a leading contribution to the integrals of the first term in (13 . 16[) in the large- A" limit. The 
saddle point equations read 



dv£\x,y) N-l 








dx h 

r(0) 



— (V{x)-y 2 -tiR x >{x)) 



(3.19) 
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where Rx(x) and Rx'{x) are the resolvents for X and X', respectively: 



In the large- iV limit, these two become coincident, and the explicit form is given by 
eqs. (IA.40I) and ( 1A.41I) in Appendix |A] Since it is seen from (1A.44I) that the origin y = 
belongs to Iy (which is also suggested by the Z 2 symmetry), the second equation in (13.191) 
leads to the solution x — (1 + c)/2 = as a saddle point of x. This value coincides 
with the critical point of x, at which a cubic equation satisfied by the universal part of 
the resolvent Rx(x) becomes triplv degenerate as Rx(xo) 3 = when g = g*, c = c*. (For 
details, see ref. [8] or Appendix lAln . The nonuniversal part of Rx(x) is given by 

RT(x) = ^ 2V '( X ) ~ + - R ^ x ) = R T{x) + Rx(x). (3.21) 

Making use of f)3.2ip to the first equation in (13.191) determines saddle points of y as 

y = ±2c = ±Q„ (3.22) 

which respects the Z 2 symmetry. From eqs. (IA.42I) . (IA.44[) in Appendix |XJ we recognize 
these saddle points to be outside the supports of the eigenvalue distributions in the double 
scaling limit □ 

x Q = P*£ I x , yo = ±Q* i Iy- (3.23) 
At the saddle points, v}h(x, y) in (13 . 1 8[) takes the form 

N — 1 / v ' 

Kff(P„±Q.) = — jj-V(P*) - (trlog(P* -X')) ■ (3.24) 

3.2.2 The second term in ( 13TT61) 

For the case y G Iy, we have the solution y Q = G Iy in the second equation of (13.191) . The 
solution is also expected from the Z 2 symmetry. However, as noticed above, in this case we 
cannot trust the form of V^\x, y) given in (I3.18p . This kind of saddle point first appears 
in the dual model, not seen in the case of the original model [HJ, [T] . In order to resolve this 
problem, we consider performing the F-integration first in the partition function and then 
calculating instanton effects from isolated eigenvalues in the resulting one-matrix model 
for X. In evaluating the second term of (I3.16p . we replace the integration over the interval 



7 Hereafter, in taking the large- N limit, c is fixed at the critical value 



_ -1+2V7 



27 



8 It turns out that the o(a) term of 7 in (|A.44[) is positive, by considering next-to-leading contributions 
to the solution of (IA.34j) . Thus, strictly speaking, y is slightly inside the cut Iy by the order of o(a), 
and coincides with the right or left edge of Iy in the double scaling limit. Hence there is some subtlety in 
justification of V^(x,y) given in (|3.18p . However, since the quantity o(a) is in fact negligible compared 
to contributions of 0(a) usually playing relevant roles in the double scaling limit, we may assume that 
the eigenvalue distribution at yo is almost zero, and that the operator trlog(j/o — Y) is not so singular 
that it can invalidate (|3.18p . 
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Jy with that over the whole real axis. Then, saddle points in the y ^ Jy region, which are 
nothing but the ones considered in (I3.23p . would give rise to an error in this replacement. 
However, as shown in Section 13.3.11 they only give contributions exponentially small by 
the factor ( 13.401) compared to those from the integration over y G ly. Thus, concerning 
the leading contribution of the second term in (13.161) . we can neglect such an error and 
justify the replacement of the integration region. 

We go back to the expression (13.61) to rewrite the second term in (I3.16P as 

N j dxdy e ~ v ^ y) 



ZN-i(h') J X 3 
Dn-i 



J \t\ Jit 'lW(W + A*i) 

dxD N [ dXdYtr8{x-X)e~% tlSD ( x ' Y \ (3.25) 



where we have inserted 1 = J d\]y5(x — \n) and used the fact that the integrand is 
symmetric under the interchange among {Aj,x} {i = 1, • • • , N — 1). The measures dX, 
dY are normalized as (13. 2J) . As the result of the Y- integration, we obtain 

N I dxdye- y ^ y) 



Z N -x(h') J X X 
x exp 



dx / dXtx5(x -X) 



jtrV(X) - itrlog (l <g> 1 - -^(X 



1 ®X) 



(3.26) 



Next, integration over the angular variables of X yields 

N j dxdy e - v ^ x ' y) 

J x<£lx,y£lY 



N 



ZN-x{h') J x £I x 



N 



N 



dxJ^ 1 / I[ rfA A (Ar) (A) 2 ^5(x-A,) 



A=l 



x exp 



i=X 

N oo / x „n N N 

i=1 n=X x ' k=0 i=l 7=1 



n—k 



(3.27) 



where the normalization constant J§ arises upon the angular integration and satisfies for 
an arbitrary (iV) -invariant function of X: f(trX, trX 2 , • ■ ■ ) 



/ dXf(tiX,trX z ,---) 
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j£ is calculated in Appendix [B] as f]B.24|) . We rename one of Aj's as x by integrating the 
delta functions Yli=i ~ ^i) an d introduce the (N — 1) x (N — 1) Hermitian matrix X' 
again to obtain the effective potential V e s(x): 



N 



dxdy e 



-V e s(x,y) 



N- 



' N 



Zn-IV 1 ') Jxpx 



dX' det(x-X 



l\2 



x exp 



N — 1 



N - 1 



h 



— trV{X') - —f^-V(x) - - log I 1 - — 1 - tr log I 1 



2x 



x + X' 
1 + c 



-trlog (10 1- ~^—{X' <g> 1 + 1 <g> X') ] 
2 \ 1+c / 



iV 



n- 1 n 7^ _1 7* 
J^N-l^N^N J N-1 



Z N -i(h') 



x exp 



x exp 



/ dx [ dX'dY' det(x - X'f (c + 
N - 1 1 

— 77— V(a;) - - log(c + 1 - 2x) - tr log(c + 1 - x - X') 
N- 1 



/r 



— tr V(X' 



1 + c, 



Y" - X'Y 



.Dm J? 



N-l 



N- 



D N -iJ? 



D N J X 



N 



dx 



detfx - X 



A2 



D 



^(o+l)"-* / 



,^ /x \det(c+l-x-X0 

rfxe -V off (x)_ 



^iy(x)-ilog( C +l-2x-) 



iV-lVjV 

Therefore, when y G Jy, the effective potential for x Ix is expressed as 

e -v B «(x) _ I det(x-X') 2 \ n^l v(x) _i log(c+ i„2x) 
\det(c+ 1 - x -X')l 

In the large-X limit, the leading term of V e g(x) is reduced to 



N-l 



y^V(x) - 2Re (tr log(x - X'))' d + (tr log(c + 1 - x - X')) 



h 



(3.29) 



(3.30) 



(3.31) 



Note that for large x ^ Ix, the third term may yield the imaginary part which prevents 
us from interpreting V^\x) as the effective potential. However, at least near the saddle 
point for x to be determined in Section 13.3.21 the point c + 1 — x is also outside the cut 
and the third term stays real. Since v£g\x) is again proportional to X — 1, we can apply 
the saddle point method to compute the second term of (13.161) in the large-X limit. 
The saddle point equation from V^\x) reads 



dx 



-(N-l) (2ReR x >{x) + R x >(c+ 1 -x)) 



(3.32) 
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where we have used (13.211) . The saddle point equation (13.321) is solved in the double 
scaling limit in the next subsection. 

Finally we comment on the third term in ( 13. 16ft . The solutions ( 13.231) for the equations 
(13. 191) are unique in the case y ^ Iy, and thus ( 13.181) has no saddle point in the region 
x G I x , y 4- Iy- Although strictly speaking the expression of the effective potential in the 
large- N limit given in (I3.18P can not be sufficiently trusted when x G I x , V ^ Iy-, here we 
assume that the third term in (I3.16P does not contribute to the numerator in (13.131) . 



3.3 Nonperturbative effects 

The leading term of the denominator in (I3.13P is computed as 



dxdy e 



-V e g(x,y) 



exp 



leading 



FP AT 1 

(trlog(/3 - X% + (N-1) R x ,{x')dx' - —r^V((3) 

Jc+1-0 n 



(3.33) 



where j3 is the right edge of the eigenvalue distribution of X in the one-matrix model 
appearing in (13.261) : 



ZnQ 1 ) = J dX exp 



~trV(X) - hr\og ( 1 ® 1 - t^—(X <8> 1 + 1 ® X) 



. (3.34) 



Since this model is obtained simply by performing the Gaussian integration over Y in the 
dual model (13.11) . the partition function is exactly the same as in ( 13.11) . Details in the 
derivation of (I3.33|) are given in Appendix [B] (See eq. (IB. 191) ). 



3.3.1 The case x G" Ix, yo & Iy 



At the saddle points (I3.23p . the leading term of the numerator in (13. 13j) is given by the 
effective potential (13.241) . Putting it together with the contribution from the denominator 
( I3.33p . the leading term of the chemical potential fi takes the form as 



(/i leading) = exp 



N - 1 



+ {N- 1) ( I R X '{x')dx 



Rx'{x')dx' 



C+1-/3 



(3.35) 



Note that the potential terms above exactly cancel the nonuniversal parts of the resolvents, 
by using the expression (13.211) . to give 



(n leading) = exp 



{N / R x >{x')dx' 



R X '(x')dx' 



C+1-/3 



(3.36) 
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It can be written solely in terms of the universal part of the resolvent, and the nonuniversal 
part does not contribute. This observation would be relevant if we consider proving the 
universality of nonperturbative effects in the dual theory as done for the standard c < 1 
noncritical string theories [H EJ [7j. 

In order to take the double scaling limit, let us introduce the scaling variable ( as 
x' — + a() with the lattice constant a. Then, the resolvent is expressed as 

4 4 

R x (x') = ataw(Q + O (at) , w(() = (( + \J?-f^ ' + U ~ \/c 2 -t) ' , (3.37) 

where a is a numerical constant and T is proportional to the cosmological constant t 
defined by 



A 4 

Ot 2 5 

2i ■ 5c 



5 = 9* (1 " a 2 *) =g*U- fl2 ^) • ( 3 - 38 ) 



s is an irrational number s = 1 + \/7, and = V5c 3 is the critical point of g. (For more 
details in the double scaling limit of the dual model, see Appendix EJ) We carry out the 
integrations in ( 13.361) by moving to the variable ( as 

x' = A <£> c = o, 

x' = /? = A (l-aV 7 ?) ^ C = - V / f, 

a/ = c + 1 - p = 2P* - p (= Vf, (3.39) 
and end up with the result 

(// leading) = exp ^-5^^y^(iV - l)aM j . (3.40) 

Since the sphere free energy is expressed as 

N 2 F = ^5^N 2 a^t 7 3 (3.41) 

in the parametrization used here^l, we take the double scaling limit in such a way that 
f)3.4ip agrees with the first term in (I2.14p . Explicitly, 

1 1 7 

iV -> oo and a -> with — = 2 • 5eiVa3 fixed. (3.43) 



9 Eq. (|3.4ip is obtained, for example, by integrating the leading term of W3 in (|A.15[) with respect to 
t, since 

in the original model, and the spherical free energy takes the same form for both of the original and dual 
models. 
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Then, ( I3.40p becomes 



(// leading) = exp ( -^-^ ) , (3.44) 



which exactly reproduces one of the leading terms in the nonperturbative effects given 
in ( 12. 161) . Moreover, we have found two degenerate saddle points ( 13.231) which give this 
nonperturbative effect. Notice that the original c = 1/2 theory also has two degenerate 
nonperturbative effects exactly given by ( 13.441) [3j [7] . Thus, the leading nonperturbative 
effects ( 13.441) precisely match between the original theory and the dual one including their 
multiplicity. 

3.3.2 The case x g I x , Vo G Iy 

From eqs. ( 13.311) and ( 13.331) . the leading contribution to /i, which comes from a saddle 
point x = xq Ix satisfying ( 13.321) . takes the form 



(/i leading) = exp 



2 J Rx>(x')dx' + J R X '(c+ 1 -x')dx'J . (3.45) 



It is again expressed only by the universal part: 



(fj, leading) = exp 



{N-l) J {2R x >{x , ) + Rx>{c+l-x'))dx' 



(3.46) 



as anticipated from (I3.32[) . Note that both of Rx>(x') and Rx>(c + 1 — x') do not have 
the imaginary parts for (3 < x' < Xq. In terms of the scaling variable (, the saddle point 
equation (13.321) is written as 

2w(() + w(-()=0, (3.47) 
and it is easy to see that a solution to this equation on the first (physical) sheet is given 



by Co = yT/2. Substituting this value into (I3.46p . we find another nonperturbative effect 

(// leading) = exp ( -5^^-y^(iV - l)^ a lj . (3.48) 
In the double scaling limit (I3.43p . this becomes 

(/i leading) = exp ^-^^-t^j , (3.49) 

to reproduce the leading term of the remaining nonperturbative effect in ( 12 . 16D . It is also 
exactly the same as the remaining one of three nonperturbative effects in the original 
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c = 1/2 string theory. Combining the results in (13.44j) and (13.491) , we conclude that the 
T-duality transformation does not change the number of the ZZ-branes and their actions 
(weights). 

Before closing this section, we make the following one comment. The disk amplitude 
(trlog(y — Y')) d (or Ry(y)) represents a quite singular dual spin configuration along the 
boundary, and its counterpart in the original model does not appear in the instanton 
calculus [7]. It is worthy of noticing that such amplitudes cancel out due to the Z 2 
symmetry in the process of the computation of IQj (x, y) in (13.171) and they do not affect 
the interpretation of the dual ZZ-branes. 



4 Coefficients of the nonperturbative effects 

In this section, we consider next-to-leading terms in the chemical potential of the instan- 
ton, namely, coefficients in front of the nonperturbative effects (13.441) and ( 13.491) . 

The coefficient coming from the denominator in (13.13!) is evaluated in Appendix [B] In 
eq. ( IB.40p . the first term in the exponential 2(iV — 1)R represents nothing but the leading 
term (I3.33p . and the remaining parts 

(2n)%y/(N-l)h' exp (r + h'^j (4.1) 

contribute to the coefficient. R is defined by eq. (1B.38j) . and R + hi '|0 can b e written in 
terms of several disk and cylinder amplitudes as (1B.43h . 

In the following, let us consider contributions to the coefficient from the numerator in 
each of the cases x$ £ Ix, yo & ly and xq Ix, Vo G ly- 

In the first case xo ^ Ix, Vo ^ Iy, the next-to-leading contributions are provided by 
the O(N ) part in the exponent in (I3.15P evaluated at the saddle point and also by the 
Gaussian integration of (I3.17P around the saddle point. The latter is easy to calculate to 
yield for both of (x ,y ) = (A, ±Q*) = ( 1 ± £ , ±2c) 

N [ dxdy e ~ v ^ y) = i— e - y c ( ffWo) e -v; ( ff Wo) (\ + o ( — J J , (4.2) 



2c V VA^ 



where 



V£\x,y) = ^(v(x) + ^y 2 -xy^-(tT\og(x-X')y d , 

V$lx,y) = -l((tr\og(x-X'))X-((tr\og(y-Y'))% 

+ (tr\og(y-Y')tr\og(y + Y% (4.3) 

are obtained from (13.151) by using the Z 2 symmetry. For both of the two saddle points (14.21) 
takes the same value, which should be expected from the Z 2 symmetry. Together with 
the contribution from the denominator, e~ v e« ( x o,yo)-2{N~i)R gj ves leading part (13.441) . 
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and Vgg(x,y) represents the next-to-leading term in V e g(x, y). The cylinder amplitudes 
appearing in (x, y) will be found by integrating with respect to z and z' 

tr — -^rrtr — -\ , / tr — \— tr — -\ , / tr — 3— tr ^— rr\ . (4.4) 

z-X' zf-X'/e \ z-Y' z' -Y' I c ' \ z-Y' -z'-Y'/ c v ; 

For this purpose, we will need the explicit form of these cylinder amplitudes without 
taking the double scaling limit. In fact, as long as the first one in f |4.4j) is concerned, 
such an amplitude is calculated in [33]. However, the result has a complicated form to 
makes it difficult to perform the z, ^'-integration explicitly. Moreover, to the best of our 
knowledge, the cylinder amplitudes for Y' like the remaining two in ( 14.4ft have not been 
found in the literature. 

Similarly, in the second case x Ix, yo £ Wi from (13.291) the numerator takes the 
form of 

N I dxe~ v ^ {x) 



(4.5) 



where x = P* (1 + «Co) = P* U + a\/ r / 2 



Kf?^) = ^ 7 ^^(a;)-2Re(trlog(a;-X / )) , d + (trlog( C +l-x-X / )) , ci , 

K?^) = -J((trlog(x-I') 2 ) 2 };^((trlog(c + l-x-I')) 2 ); 

+ (tr log(x - X') 2 tr log(c + 1 - x - X')}' , (4.6) 



and e 0"°) 2 ( N 1 ) R provides the leading contribution (13.491) . The factor coming from 
the Gaussian integration can be evaluated easily in the double scaling limit. Then we find 



N I dxdye- y ^ = J^e-^We^'W. (4.7) 

2 

Note that the factor has the a-dependence of a~3 and it is real. Therefore, it is important 
to check that e~ VaS <x °- 1 exactly cancels the a-dependence in the prefactor so that the 
nonperturbative effect will be finite in the double scaling limit. Also, it is interesting to 
see whether e _v °ff ^ provides the factor i. In examining V^\xq) } however we encounter 
again the technical problem in the z, ^'-integration of the first amplitude in (I4.4p . 
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5 Conclusions and discussions 



Based on the dual two-matrix model, we formulated the chemical potential of the instanton 
and showed that the instantons in the dual model give the same nonperturbative effects as 
those in the original model in the leading order. Since the dual model defines the T-dual 
theory of the standard c = 1/2 noncritical string theory, it implies that we have identified 
the degrees of freedom of the dual ZZ-branes, which provide nonperturbative effects in 
the dual theory, and shown that the number of species of the dual ZZ-branes and their 
tensions coincide with the original ones. Thus our result gives further support on the 
identification of the ZZ-brane or nonperturbative effect in string theory as the instanton 
in the matrix model. 

As mentioned in Section [21 this result is expected because the free energy of the dual 
model is the same as the original one (rigorously speaking, in the sense of the genus 
expansion) and hence it should satisfy the same string equation in the double scaling 
limit as the free energy of the original theory does, which then implies that both theory 
have the same form of the leading nonperturbative effects. Another argument which 
supports our result is that the leading part of the nonperturbative effect is given by the 
disk amplitude in the ZZ-brane background, but at the level of the disk topology the 
T-dual relation (12.91) holds. In this interpretation, it is crucial that the disk amplitudes 
(trlog(y ± Y'))' d cancel out in the process of the computation and do not contribute to 
the instanton effects as shown in Section [3J 

However, regarding the next-to-leading part in the chemical potential of the instanton, 
the string equation does not guarantee that the dual theory has the same coefficients in the 
nonperturbative effects as those in the original theory. In the direct calculation from the 
matrix model in Section HI several cylinder amplitudes appear to contribute to them. In 
particular, the amplitudes ((trlog(y — Y')) 2 )' c and (trlog(?/ — y')tr log(y + Y'))' represent 
quite singular dual spin configurations and their counterparts in the original model do not 
appear in the calculation of the nonperturbative effects [7J. Thus, if the above amplitudes 
persist in contributing after taking the double scaling limit, there is a possibility that 
the T-duality does not hold at the level of the next-to-leading order. From this point of 
view, the comparison of the coefficients in the nonperturbative effects in both theories is 
quite interesting because it may possibly reflect the violation of the T-duality. On the 
other hand, if the nonperturbative effects computed from the original and dual matrix 
models coincide even including their coefficients, this may imply that the matrix model is 
more fundamental than the string equation and that the nonperturbative effect in string 
theory has large universality which contains the T-duality. In fact, it is shown in [7] 
that the coefficient in the nonperturbative effect is universal within the original c < 1 
string theory in the sense that it does not depend on details of the potential of the matrix 
model. Therefore, it is an interesting problem to investigate whether this universality 
involves even the T-duality. We hope that this kind of study would give some insight 
into universality of string theory itself. For this purpose, it is desirable to calculate the 
cylinder amplitudes of trace-log operators like (14.31) and (14.61) which contribute to the 
coefficients. We hope we will be able to report on it in the near future. 
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Apart from the coefficients of the nonperturbative effects, there are lots of quantities of 
interest to be computed. For example, in order to reinforce the identification between the 
dual ZZ-brane and the instanton in the dual matrix model, it is preferable to calculate the 
loop amplitudes under the corresponding backgrounds for both of the continuum theory 
and the matrix model and then to compare the obtained results. In particular, it is 
quite interesting how the T-duality and the dual ZZ-branes are realized in the continuum 
Liouville theory. In the case of the standard c = 1/2 conformal field theory, we know 
how the T-duality maps the relevant operators into themselves. If we clarify it for the 
gravitationally dressed case (c = 1/2 noncritical string theory), we will be able to make 
more explicit the correspondence between relevant operators and associated ZZ-branes in 
the dual theory and to construct the T-duality transformation law of the ZZ-branes in the 
noncritical string theory. The result in this paper strongly suggests that the T-duality 
transformation rule of the relevant operators is the same as in the nongravitational case 
and that the T-duality does not change the tension of the ZZ-brane associated with each 
relevant operator. 

Finally, in refs. [25J and [8J, c = 1/2 noncritical string field theories corresponding 
to the original and dual two-matrix models have been constructed respectively, and it is 
clarified how string fields in the two theories are related each other under the T-duality. It 
is quite interesting to express the original and dual ZZ-branes in terms of the string fields^ 
and to determine the T-duality transformation rule between the ZZ-branes by utilizing 
the established transformation rule of the string fields. It will be also intriguing from 
the viewpoint of finding the relations between different descriptions of nonperturbative 
objects in string theory - (ZZ-)branes and string fields. 
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A Double scaling limit and Resolvents in the dual 
model 

In the paper [8], various amplitudes in the dual two-matrix model (12. 7p are calculated and 
their expressions in the continuum limit are obtained. Here, for convenience, we present 
the result of the disk amplitude for the resolvent operator of X: 




(A.l) 



For somewhat related works, see refs. [33]. 
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derived in [5]. Also, we compute the disk amplitude for the resolvent of Y: 

which has not appeared in the literature. 

Since the dual model (12 .TP is obtained from the original one ( 12 .ip by changing the 
integration variables as (12. 5p . the partition functions for both models coincide and have 
the same critical point. The original model has the critical point 

j Q v// 

c* = 27 ' 9 * = V^^' ( A - 3 ) 

and the double scaling limit to the c = 1/2 string theory is taken as iV — ► oo and the 
lattice spacing a — > with fixing 

c = c*, t=— — ^ and — = (const )Na*. (A.4) 
9*a 9s 

t and g s represent the cosmological constant and the string coupling constant in the 
continuum theory. Correspondingly, the critical point of the dual model is 

(c, g) = (c, &,) = ^ 1 27^ ' ' ( A ' 5 ) 

and we obtain the dual c = 1/2 string theory by taking the double scaling limit iV — > oo, 
a — >• with keeping 

= 0*, i = ^ — and — = (const)A^a 3 (A. 6) 

9*a 2 g s 

fixed. In what follows, the coupling c is understood to be fixed at the critical point c*. 
A.l Rx(x) 

As in [8], [9] , the following closed equation for Rx(x) is derived by combining five Schwinger- 
Dyson equations obtained in the planar limit: 

gRx{xf + f 2 Rx(x) 2 + fiRx(x) + f = 0, (A.7) 



where 



h = g 2 x 2 + {5c-l)gx-2c(l + c), 

f l = Acg 2 x 3 - (6c - 2c 2 )gx 2 + 2c(l - c 2 )x - 2g 2 V x + g(l - 3c), 
f = (6c-2c 2 -4cgx)gV x -4cg 2 V X 2-4cg 2 x 2 

+ (6c-2c 2 )gx-2c(l-c 2 ) - g 2 . (A. 



20 



We used the notation: Vx^y™ = (^ftr(X n F m )) . The amplitudes Vx™ (n = 1,2) are 
related to the amplitudes in the original model W n = ^trA n ^ d (n = 1, 3) as 



After the shift 



(1A.7j) becomes 



with 



V x = V2W U V X 2 = - — —W 1 - -W 3 - -. 

eg c c 



f 2 

Rx{x) = -— + Rx(x), 



R x (x) 3 - ^Rxix) -^F = 
F\ = -7^/2 — ~fi, Fo = — ~/o- 



9 9 

At g = g#, F± and Fq can be written as 



9' 



9 



9 



(A.9) 
(A.10) 

(A.ll) 
(A.12) 



F 1 = -(z-s) 



3 

c\ 2 



, z — s \z — s 



+ 3^ ) (z-s- 3VQ 



(A.13) 



where z = Vbcx, s = 1 + \/7, and the fact was used that W\, W 3 expressed in [25J: 



Wx 



-8p^ + 3(2^) 2 , 2 -136pt + 27(20*) 2 



(A.14) 



3(2^) J 



+ a' 



27(20, 



^3 



W 3 non 



27(20,)' 



: (5t)s +a~ 



81(20*) ; 



-(5t)3+a 



-8527/4 + 972(20^ 



972(20,)^ 



■t 2 + a~ 



W™ n + W 3 , 

32(420 - 839p*)p| + q2 160(252 - 611p,)p^ 



(A.15) 



as- 



729(2^) 5 
320p^ 



729(2#*) E 



,„ ,4 10 160p° 
-(5t)3+ a - 



70(1152 -3593p*)p* 



81(20*) 5V ~" y ' " 243(2^)5 
with p* = 3c become at the critical point 



W 3 * 



9 

40c 
32/9 



40 

3c H 

27 



(140 - 839c) . 



243 V40c / 

Eqs. (1A. 13j) determine the critical point of x, denoted by P*, as 

Pi = 



5c 



(A.16) 
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where the cubic equation (lA.lip becomes triply degenerate: Rx(Pi) 3 
equation (lA.lip is solved by 



Rx(x) 



F + v/F 2 - 4Ff 



+ 



F - v/F 2 - AFf 



0. In general, the 



(A.17) 



Introducing the variables in the continuum theory T and ( as 



g = g*[l-a — T 



x 



(A.18) 



(T is proportional to the cosmological constant t), the continuum limit of the solution 
(IA.17D takes the form 



R x (x) = a*a'w(0 + Oia*), w(() = [( + VC 2 - T 



with a' 



1 4 



. Also, 



R x (x)= R^ n (x)+R x (x), R 



c-vc 2 -t 



h 



(A.19) 



(A.20) 



where i?™ n and Rx stand for the nonuniversal and universal parts, respectively. 

Eigenvalue distribution of the matrix X is seen as a cut of Rx{%) on the x-plane, 
which is determined in principle from (IA.170 under the condition of Rx{x) having one 
cut and being analytic for Rex > x c with x c some finite number. However, practically 
it is not an easy task in general, since _F 2 — 4Ff is a complicated polynomial of x with 



the degree ten. w(() in (1A.19I) has one cut £ 6 ( — oo 



, which means that the right 



edge of the cut for Rx(x) is x = (3' = PI — ay Tj when g is near the critical point g*. 

The leading value of the left edge also can be determined from the expression of (1A.13I) . 
Plugging it into flA.171) . we find the form of Rx(x), with c and g set to the critical values 
but x left generic, as 



Rx{x) 



1 

3 V 5 



(A.21) 



g±(z) = — (z — s + 3y/6j (z — s — 3\/6j ± 6y/3~J— (z — s + 
Here, it is easy to see that g+(z) has no zero and g~{z) has the unique zero at z = s. For 



z ~ s, 



g+(z) = 108 + 0((z-s) 2 ), g_(z) 



1 



108 



s) 4 + 0((z-s) 



(A.22) 



From these properties of g±(z), it can be shown that ( 1A.21|) has one cut of the interval z G 



[s — 3\/3, s]. Thus, the left edge of the cut for Rx{x) is given by x = a' = Pi — ^= + 0(a 
in the case of g near the critical point. 
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e -NtrS D (X,Y) 



e -NtrS D (X,Y) 



e -NtrS D (X,Y) 



-NtrS D (X,Y) 



A.2 R Y (y) 

To obtain a closed equation for Ry(y), we combine the six Schwinger-Dyson equations 
derived from the following identities in the planar limit: 

= / rtd^y V An f 1 ^ e-^tr to 

j a dYa \y~ Y J 

o = /ft/rE^^r) 
o = /A^EAt,^) 

= J d N2 Xd N2 Y^^tr(^— l —t a YX^e- mr ' s ^ Y \ (A.23) 

Here, (a = 1, • • • , A^ 2 ) are a basis oi N x N Hermitian matrices satisfying 

tr(Wt a Zt a ) = tr W tr Z, tr (^) tr(Zt Q ) = tr(WZ) (A.24) 

a a 

for arbitrary matrices W, Z. The result is 

R Y (y) 4 + hRyiy) 3 + b 2 R Y (y) 2 + b 3 R Y (y) + b 4 = 0, (A.25) 

where 

61 = -8cy, 

62 = 4^V+(-l + 2c+19cV + 2(l + c-2^ x ), 

6 3 = -16^V + 4c(l + c)(l-3c)y 3 -8c(l + c-2^Vx)y, 

64 = 16c#V 

+{-4c(l + c)(l - 3c) - 4(1 - c)(l - 3c)^Vx + 4(1 - 3c)^V x2 + 8^V xy2 }|/ 2 
+ (l + c-2£Vx) 2 . (A.26) 

(Note that Vy = Vxy = from the Z2 symmetry under F — > —V.) Also, Vxy 2 can be 
expressed in terms of the amplitudes W\ and in the original model, similarly to Vx 

and Vx^, as 

v/2 V. c^ 2 c c^y 
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We shift R Y (y) as 

R Y (y) = -\b 1 + R Y (y) J (A.28) 
so that the cubic term of Ry(y) in (1A.25I) vanishes, to obtain 

R Y {yf + B 1 R Y (y) 2 + B 2 R Y {y) + B 3 = 0, (A.29) 

where 

Bi = -\b\ + b 2) 

1 , 1 
B 2 = -bi - -hb 2 + 6 3 , 

o z 

S 3 = -^+^-^3 + t 4 . (A.30) 
25o lb 4 

The critical points of y are determined by 

£^ = 5^ = B 3 l = 0, (A.31) 

with the symbol "I*" meaning c, g, Wi and W3 set to the critical values. The equation 
for B 2 is trivially satisfied, and the remaining two lead to 

Bil = 20c 3 (y-g / J 2 ( 2 / + Q / J 2 = 0, 

B 3 l = -80c 5 (y-g'J 3 (y + g'J 3 = 0, (#= * V (A.32) 



. 5c 

to determine the critical points as 

y = ±&. (A.33) 

The result reflects the Z 2 symmetry as should be. We introduce the scaling variable 77 as 
V — 01(1 + ar ?) an d take the continuum limit of the quartic equation flA.290 by using the 
expressions of W\, W3 in I I A. 141) and ( 1A.15I) . The result is 

R Y (y) 4 - a 2 ^c(32 V 2 - s 2 f )R Y (y) 2 - a 3 2 0^ rf + O (a^) = 0, (A.34) 



which is solved as 



R Y (y) = EST(V) + MV), 

RT\y) = -\h = 2cy, 



R Y ( y ) = a j 8 • 23 y|r/t + O (a* ) . (A.35) 
(We took the branch of Ry{y) to be real positive for y or rj large positive.) The scaling of 

a 3 

Ry{y) with the power a 4 is clearly different from the scaling of correlators ever computed 
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in the original and dual models [HI [22] , which is characterized by the power of cubic roots. 
Ry(y) is regarded as a disk amplitude with the disorder operators distributed densely over 
the boundary, and such an amplitude has not been computed in the literature concerning 
either of the original and dual theories. Also, the leading expression of Ry(y) does not 
depend on the cosmological constant, meaning that the disk collapses into a surface of 
zero area although its boundary has a finite length. In the above, we took the continuum 
limit around the critical point y = Q'*- Similarly, the limit around the other critical point 
setting y = — Q'*(l + arf) leads to the same equation as (IA.34[) . which is a consequence 
from the Z2 symmetry. 

The expression ( 1A.35|) has one cut r\ G (—00, 0], whose end point does not depend on 
T. It shows that the right edge of the cut y = (j* on ?/-plane receives no 0(a) correction 
even when g moves slightly off the critical point. In order to see the global structure of 
the cut on y-plane, we solve eq. (IA.290 by using (IA.320 in the case of (c, g) fixed at the 
critical point but y left generic. The result is 



Rr{y) = VlOJ>(y-Qi)Hy + Q'* 



y-\^(y-Q0(y + Q'*) 



(A.36) 



(We took the same choice of the branch as above.) which has one cut of the interval 
y G [— Q'*,Q'*\- Therefore, we see that Ry(y) has the Z 2 -symmetric one cut of the form 
y G [—7', 7'] with 7' = + o(a) for the case g near the critical point. 



A. 3 Expressions for the action (13.11) 

The calculation of the instanton effects in Section [3] is done by using the rescaled action 
Sd(X, Y) in (13. ip instead of Sd(X, Y) in (12.71) . The former is obtained from the latter by 

X -> \ X, Y — > -7- Y. (A.37) 
9 9 

Here, we consider the effect of the rescaling to the results obtained in Sections IA.1I and 
IA.2[ to match them with the expressions in Section [3j 

Let us rename the resolvents calculated using the action S D (X,Y) as Rx(x) and 
Ry{y)) which are nothing but those appearing in the sections |A. II and IA.2I As a result of 
the rescaling, the resolvents Rx, Ry considered in Section [H] are related to Rx, Ry as 

Rx{x) = Ir x (Ix), Ry(y) = \k Y (\y\. (A.38) 
9 \9 J 9 \9 J 

Hence, the critical points of x and y for Rx{x) and R Y (y), denoted by P* and Q* respec- 
tively, are given as 

P* = g*K = sc= -4^, Q* = g*Q'* = 2c. (A.39) 
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For Rx{x) = R^ n (x) + Rx(x), the nonuniversal part becomes 

1 (\ \ 1 



RT(x) 



1 
— X 

9 



^(x 2 + (5c- l)x-2c(l + c)) 



1 

3/i 



^'(x) - V"(l + c-x)), 



(A.40) 



where /i = g 2 , and = — |x 3 . The scaling variable ( is introduced as x 

A(l + and the universal part takes the form 



Rx(x) 



\r x (\x 
9 \9 



a3&w(C) + Ola's), 



n 4 

S -i 



a 



9* 



-a 



2a -5c 



™(C) 



c + vc 2 -t 



c- vc 



T 



(A.41) 



When <? is near the critical point, endpoints of the cut of Rx(x): x 6 Ix = [ct,0\ are 
expressed as 



P 
a 



g*0 = P.(1- aVT^ 
g,a' = A -3v^c + C?(a). 



(A.42) 



Similarly, for Ry(u) = -Ry on (y) + Ry(y), we define the variable r] by y = Q*(l + aq), and 
obtain 



9 \9 



l Rl on (-y 



1 

4£ 



2c 



My) 



9 \9 



i 



3 8-24 3 /n . 
(24 — : ^4 + O I Q4 



5c 



The cut of Ry{d) is given by the Z2-symmetric interval y G ly as 

J y = 7 ] with 7 = + o(a), 



(A.43) 



(A.44) 



for the case g near the critical point. 

Since the one-matrix model (I3.34p is obtained after the Gaussian integration over Y 
in fl3.ll) . the partition functions of (13.1 1) and (13.341) are identical and have the same double 
scaling limit (1A.6j) . Also, concerning the correlators among operators independent of Y 
(for example, Rx(x)), both models give the same result. 
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B Computation of denominator 



In this appendix, we compute the denominator in (13.13!) based on the method developed 
in [32]. Using ( 13.121) . the denominator can be written in terms of the partition functions 

as 

D N Z^ st \h) 



L 



dxdy e 



-V eS (x,y) 



(B.l] 



The difference between the O-instanton partition function mt,t \h) and the total par- 



tition function Z^{K) is exponentially small as exp(— C/g s ) which is nothing but the 
nonperturbative effect. Since it is negligible in the computation of fi, we can replace 

^(0— inst) i 



, v \h) with Z N (h) in (JET]). Therefore, 

dxdy e- Vcs{x ' y) 



/ 

J (x 



D N Z N (h) 
D N ^Z N ^(h') 



(B.2) 



'(x,y)&S 

Namely, the denominator can be obtained basically as the ratio between the partition 
functions of the matrix model with rank N and that with rank N — 1. As mentioned in 
( 13. 3p . under the measure (I3.2p the total free energy has the standard 1/iV-expansion 



Z N (h) 
Z N -!{h') 



exp 
exp 



1 



-N 2 F (h)-F 1 (h)+O[ — 



.(N-lfF (h')-F 1 (h')+O 



(N-iy 



(B.3) 



Hence, 



dxdy e 



-V eB (x,y) 



{x,y)eS 



D 



N 



X 



exp 



{N - l)(2F (/i / ) + tiF^ti)) - ( F (h') + 2h'F^ti) + ^h' 2 F'^h" 



1 + {n 

The computation is essentially reduced to finding the sphere free energy F (h). 



(B.4) 



B.l Sphere free energy 

In order to compute the sphere free energy, we first perform the integration with respect 
to Y in (13. ip and express Z^iK) as the integration over the eigenvalues of X. Then we 
obtain 



Z N {h) 



jx-i 



N 



exp 



N 



5>g(A ! ,-A J ) 2 --$>(A 



i<j 



55> 



AT 2 

{c + l-Xi- Xj) + — log(c + 1) 



(B.5) 



27 



where is defined in ( I3.28P whose explicit form is given in Appendix IB.4.11 In the 
large-iV limit, by introducing the eigenvalue distribution p(A) = (j^'^2 i 5(X — A*)), this 
becomes 



Zf{h) = J§ exp 



d\p{\)V{\) 

. (B.6) 



iV 2 <j J dXdpp(X)p(p) log | A — fi\ — ^ y 

l y 1 

-- / dXdpp{X)p{p) log (c + 1 - A - p) + - log(c + 1) 

Therefore, the sphere free energy is expressed in terms of the eigenvalue distribution as 
F (h) = - J dXdpp{X)p{p) log |A - p\ + ~ J d\p(\)V(\) 

+ l - J d\dpp{\)p{p)\og{c+l-\-p)- l -\og{c+l) + ^\ogjX. (B.7) 

B.2 Identity of the resolvent 

Here we derive an identity of the resolvent Rx{z) introduced in (13.201) which plays an 
important role in the following. In terms of the eigenvalue distribution p(X), the resolvent 
is given by 

R x {z) = J d\^. (B.8) 

Then the support of the eigenvalue distribution is characterized as the cut of Rx(z) on 
the complex z-plane. 

On the other hand, in the large-iV limit p(A) satisfies a saddle point equation which 
follows from (IB. 71) 

= 2P [ dp^- - \v\X) + [ dp- M 



A — p h J c + 1 — (A + p) 

= Rx(X + iO) + Rx(X-iO) + Rx(c+l-X)-^V'(X), for A on the cut. 

h 

(B.9) 

Note that, for A on the cut, A < P* as shown in ( 1A.42I) and c + 1 — A > P* is always 
outside the cut. Hence the third term in the last equation in ( IB. 91) does not have the 
imaginary part. A particular solution to this equation is provided by 

RT(z) = ^V'(z) - Hi + c-z)). (B.10) 
Therefore, the resolvent Rx(z) is in general given by 

R x (z) = RT(z) + Rx(z), (B.ll) 
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where Rx{z) satisfies 

= Rx{z + iO) + R x {z-iO) + R x {l + c-z), for z on the cut. (B.12) 



In fact, this form of the resolvent coincides with the result in I I A. 401) and ( 1A.41I) in Ap- 
pendix [A] 

B.3 Fo(h) in terms of the resolvent 



Now let us express each term in ( IB. 71) in terms of the resolvent. From ( IB. 81) . it is easy to 
see that the following identities hold: 



J dzdz p{z)p{z) log \z — z'\ 
= / R x {z)dz +logA+ / dzp(z) / ReR 

J A J a J 



'z')dz' 



dzdz' p(z)p(z') log(c + 1 — z — z) 



R x {z')dz' + logA 



dzp(z) 



c+l-z 



Rx{z')dz' 



(B.13) 



where A — > oo limit is understood to be taken eventually, and (3 is the right edge of the 
cut explicitly given in (lA.42j) . Substituting these into (IB.7[) and using the saddle point 
equation (IB. 91) leads to 

F (/i) = -~/ltrlog(/3-X)\ ~\J Rx ^ dz ' + ^h V ^ + ^h J dz P( z ) V ( z )+<?> 

(B.14) 

where J represents terms independent of the potential V 



J=-~\og(c + l) + ^\ogjX. 



(B.15) 



According to (1B.4I) . we also need derivatives of Fo(h) in computing the denominator. In 
order to make explicit /i-dependence of the sphere free energy, we rewrite F (h) again in 
terms of the integration with respect to the eigenvalues as 



F (h) - J 
1 



log 



N 



i=l 



l + C 



ifl N 

2 



1 ' 2 n( i - T ^^+^) 

i,j=l V 



(B.16) 



where we have used (IB. 5|) . This shows 
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(B.17) 



From flRl4|) and (jETTjl . for the (JV- 1) X (JV- 1) matrix model with ft' defined as flU 
we obtain 

2F (ti) + tiF^(ti) 



(B.18) 



where J7"' is obtained by replacing N and ft by iV — 1 and ft,' in (]B,15j) . We will see below 
that the last term in this equation together with the overall factor D^/Dm-i in ( IB. 41) 
gives subleading contributions not affecting the leading term. Therefore, we find that the 
leading term in the denominator is given as 



dxdy e 



-V e s(x,y) 



{x,y)eS 

exp 



leading 



(trlog(/3 - X% + (N - 1) / R x ,(x')dx 

Jc+1-0 



ft' 



(B.19) 



B.4 Next-to-leading contributions in the denominator 

In this subsection, we evaluate next-to-leading contributions in the denominator in ( 13. 13ft . 

B.4.1 Derivation of Dn and J* 



First let us evaluate the measure factors D N) J§ introduced in (13.41) and (13.281) . 

In order to calculate J§, we make a connection between the measure dX defined in 
( 13. 2D and the standard measure 



dX = Yl dX, d Y[ 2(d ReXij) (d ImXy). 



^B.20) 



i<j 



For this purpose, it is sufficient to consider the Gaussian integration as follows. In the 
case of the standard measure, 



dXe -^ x2 = (2 



7T 2 . 



(B.21) 



Comparing this to (13.21) yields 



dX 



[l-c)N 
2vrft 



dX. 



(B.22) 



On the other hand, by using the method of orthogonal polynomials, 



. N-l N 

/ HdX, ^(AfriS,^ = N\ J] K = (27r)f Y[ p \, (B.23) 

J i n=0 p=0 
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where we have used the well-known fact that in the Gaussian case the orthogonal poly- 
nomial is nothing but the Hermite polynomial with h n = \/27m!. Thus we find 

(N \ — 
£') (<^) ' • (R24) 

Similarly to J§ in (13.281) . if we introduce Jjj by 



with dY defined in (13.21) . we obtain 



JV 



VP=0 



Using these, (13.11) becomes 

dXdYe-X v ^ +1 -¥ Y2 - XY2 ) 



(B.27) 



where 



/(A, /i) = y df/exp ^trpfC/Y^J . (B.28) 



It is calculated by the method in [2H| as 



iV(iV-l) /Ar _ 1 v jv , 

iV\ 2 f W A det,-,- eT^j 



Substituting this into (1B.27[) . we find 

iV(iV-l) / N 

-<■**)-"© 1 (n 
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Comparing this equation with the definition of ( 13. 4p . we finally obtain 

N(N-l) / \r \ — 1 / at \ JV JV 2 



Thus the overall factor in ( IB. 411 becomes 

(^1^^-^(1 + 0(1)), (B.32) 



which agrees with the result in [7j. 
B.4.2 Other contributions in (IB. 41) 

Using the Euler-Maclaurin formula, we evaluate the following quantity 



log ( f[p\] = (n*" +1 " fc ) = {J1 ^^ g {N+l) -~ A N>-^ + 0(1), (B.33) 

p=0 / \fe=l 



2 ov ' 4 2 
which is also derived in [3~o]. From (1B.24D . we see that J given in (IB. 15[) becomes 



c)N J \N 



Therefore, the last term in ( IB. 181) gives 



Next, let us examine the O(N ) contribut ion in the exponent in (IB.4[) . From (IB.17[) we 
have 

h'F^(h') =-±J dzp(z)V(z) + \. (B.36) 
Combining this equation with (IB. 181) and (IB. 351) . we find 

F (h') + 2h'FW) + ±h*F"(h') = -R- ti-^R - \ log i^f) , (B.37) 
where R is given as 
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Also, (IB. 18[) is written as 

2F (h') + h'F^ti) = -2R - log (j^-j ~ L (B.39) 

We substitute flB~32l . flB~39l and flB~37l) into flB~4j) . to obtain a simple expression of the 
denominator in terms of R: 

[ dxdye-v^'ti = (2?r) VOW - l)h' exp ( 2(JV - 1)R+ R + h'^) . (B.40) 

The first term in the exponential 2(iV — 1)R represents the leading part (IB. 191) . and the 
remaining gives the next-to-leading contributions. A similar expression is obtained for the 
denominator in the definition of the chemical potential of an instanton in the standard 
c < 1 noncritical string theory [7j. 

Finally, we have a comment on the fact that O(N ) part in the exponent in the 
denominator (1B.40|) can be written in terms of a cylinder amplitude. Noting that 

h'jfi (O)' = ^ (OtrV(X% onn , (B.41) 

where the subscript "conn" represents taking the connected part of the correlator, we find 
that 

fin ii i 

h' W = y{^og{P-X>)tvV{X%- — (tr log(c + 1 - /3 - A')tr\/(X'))' C + — V((3) 

1 1 V 1 / 1 1 V 1 ,,dP 

—l tT J^c) d + 2 \—l tr c + l-[3-X' ) d ~ 2hf V ) h 8h» 

(B.42) 

where the last term vanishes due to ( IB. 91) . Thus 

R+ "'§ = <F^i trlo ^- x ''>rKi^T trlog(c+1 " /5 - x ' ) >l 

+y (tr logOS - X')txV(X% - ^ (trlog(c + 1-/3- X')tiV(X% 

+ logA- -log A'. (B.43) 

Namely, once we calculate the cylinder amplitude (to z l X ' ^ r z'-x' ) ' we snom d determine 
the O(N ) coefficient in the denominator. In fact, the explicit form of this cylinder 



+ 
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amplitude is given in [33] for the 0(n) model with arbitrary n. (For the case |n| < 2, see 
section 3.4 in the first paper of [23] •) Using that expression in our case n = 1, we find 
that the O(N ) part in the exponent takes rather a simple form 

R + ft '§7 = C dz 7S G{z) + log A " \ H dz 7S G{z} ~ \ loe A ' (R44) 

h' h f 

where G(z) is a function introduced in [33], which is universal in the sense that it is 
uniquely determined only by the homogeneous saddle point equation (IB. 121) and by spec- 
ifying its behavior near the edge of the cut and at the infinity. For details in the case 
\n\ < 2, see sections 3.2 and 3.3 in the first paper of [53] . It is quite interesting that, even 
in the case of c = string theory defined by the one-matrix model, the denominator in 
the chemical potential of an instanton also has the next-to-leading term given by G(z) 
for the O(0) model as 

dR f 13 1 3 - a 

R + h'— = / dz -^G(z)+logA = log^— - -. (B.45) 
dh> J A 4 

It would be intriguing to examine whether this property holds for other 0(n) models. In 
our case (the 0(1) model), however complex form of G(z) prevents us from performing the 
^-integration explicitly in (IB.44|) and we have not yet succeeded in obtaining a concrete 



value of the coefficient in the denominator. It would be natural to expect that it is 
somehow related to the length of the cut as in the case of the 0(0) model ( 1B.45I) . 
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